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A method to efficiently compute, in a automatic way, helicity ampli- 
tudes for arbitrary scattering processes at leading order in the Standard 
Model is presented. The scattering amplitude is evaluated recursively 
through a set of Dyson-Schwinger equations. The computational cost of this 
algorithm grows asymptotically as 3", where n is the number of particles 
involved in the process, compared to n\ in the traditional Feynman graphs 
approach. Unitary gauge is used and mass effects are available as well. Ad- 
ditionally, the color and helicity structures are appropriately transformed 
so the usual summation is replaced by Monte Carlo techniques. Some re- 
sults related to the production of vector bosons and the Higgs boson in 
association with jets are also presented. 

1. Introduction 

Multi-particle and multi-jet final states are of great importance at the 
TeVatron and at the future LHC or e^e~ Linear Collider. They serve both 
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Fig. 1. An example of common parts of the amplitudes (in red) for the e+e — > 

e"'"e~e'^e" process. 

as signals and as important backgrounds to many new and already discov- 
ered physics channels. As an example the production and decay of top 
quarks, Higgs boson(s) or SUSY particles can be mentioned. A typical 
background is the production of weak vector bosons in association with 
jets. Among others the proper evaluation of the eight jet QCD background 
will be needed. To describe the production process of a number of particles 
the corresponding amplitudes have to be constructed. This usually results 
in a very large number of terms, such that their automated construction 
and evaluation becomes the only solution. Apart from handling the number 
of amplitudes, which grows factorially with the number of particles, the in- 
tegration over the multidimensional phase space of the final state particles 
represents a formidable task. In the past years various solutions to deal with 
these problems, implemented as different codes, have been presented. Either 
they are based on traditional methods of constructing Feynman diagrams or 
alternative methods with recursive equations are implemented [1^13]. The 
new formalism based on the Dyson-Schwinger equations recursively defines 
one-particle off-shell Green function. It does not involve any calculation of 
individual diagrams but various off-shell subamplitudes are regroupped in 
such a way that as little of the computation as possible is repeated. On the 
contrary, in the traditional approach, the same parts of different Feynman 
diagrams are recalculated all over again, see FigQ increasing the number of 
steps that should be done in order to get the full amplitude. The recursive 
approach significantly decreases the factorial growth of the number of terms 
to be calculated with the number of particles down to 4" or 3*^ ^. 

Some examples of automatic parton level generators for any processes in 
the Standard Model are e.g. CompHEP [14-16], MadGraph/MadEvent [17,18], 
AMEGIC++ [19] and the HELAC/PHEGAS package [20,21]. Codes designed for 



^ To reduce the computational complexity down to an asymptotic 3", each four-boson 
vertex must be replaced with a three-boson vertex e.g. by introducing an auxiliary 
field represented by the antisymmetric tensor iif^", see [13] for details. 
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specific processes are e.g. GRACE [22,23] as well as Alpgen [24]. Very recently 
also on shell recursive equations have been proposed [25,26]. However, event 
generators based on this new method are not publicly available yet. 

In this article the algorithm based on Dyson-Schwinger recursive equa- 
tions is briefly reviewed. It has been implemented as a new version of 
the multipurpose Monte Carlo generator HELAC in order to efficiently ob- 
tain cross sections for arbitrary multi-particle and multi-jet processes in the 
Standard Model. 

2. Dyson-Schwinger Recursive Equations 

Dyson-Schwinger equations give recursively the n— point Green's func- 
tions in terms of the 1—, 2—,. . ., {n — 1)— point functions. These equations 
hold all the information for the fields and their interactions for any number 
of external legs and to all orders in perturbation theory. The recursive con- 
tent of the Dyson-Schwinger equations for QCD has already been introduced 
in Ref. [13] and reviewed recently in Ref. [27]. To include the electroweak 
sector, new vertices for leptons, the vector gauge bosons as well as for the 
scalar Higgs boson must be included. Additionally, the recursive equation 
for (anti)quarks should be rewritten to express their interaction with the 
electroweak gauge bosons. 




Fig. 2. Recursion equations for the Higgs boson. 

In order to better illustrate this idea let us present as an example the 
recursive equations for the Higgs boson interaction with massive particles 
only. Let pi,p2, ■ ■ ■ ,Pn represent the external momenta involved in the 
scattering processs taken to be incoming. For a vector field we define a 
four vector bJ^{P), which describes any sub-amplitudes from which a vector 
boson V with momentum P can be constructed. The momentum P is 
given as a sum of external particles momenta. Accordingly we define a 
four-dimensional spinor ip^[P), which describes any sub-amplitude from 
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which a fermion with momentum P can be constructed and by tp^iP) a 
four-dimensional antispinor. Additionaly we have to introduces a scalar 
H{P) for a Higgs boson. The content of Dyson- Schwinger equations in this 
case can be understood diagrammatically as in Fig|2 The subamplitude 
with an off-shell Higgs boson momentum P has contributions from three- 
bosons and four-bosons vertices plus the fermion antifermion vertex. The 
black blobs denote subamplitudes with the same structure. The following 
general recursive equations can be written down for a Higgs boson with the 
momentum P: 

H{P) =j2KP-Vi)H{pi) 
1=1 

+ ^ShvyHh 6;(pi)6^'^(p2)e(pi,P2) 

P=pi+P2 
P=Pl+P2+P3 

+ Y ^Shff Hh 1p''{piW{p2)e{pi,P2) 

P=Pl+P2 

+ Y ^5hhh Hh H{pi)H{p2)e{pi,P2) 

P=Pl+P2 

+ Y ^5hhhh Hh H{pi)H{p2)H{ps)e{pi,P2,P3) 

P=Pl+P2+P3 

where the Higgs boson propagator is given by 

Hh = -2 2^—^ (2-1) 

P'^ — — iL h"^h 

and e{pi,P2,P3) = il is a sign function, which takes into account the sign 
change when two identical fermions are interchanged. 

The scattering amplitude can be calculated by any of the following re- 
lations, 

bJ^{Pi)b^^{pi) — vector bosons 
H{Pi)H{pi) — Higgs boson 



A{pi, . . . ,pn) = < 



(2.2) 



ip^ {Pi)ip^ (pi) — incoming fermion 
ijj^ {pi)ijj^ {Pi) — outgoing fermion 
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where 

Pi = Y^Pj^ 

so that Pi+pi = 0. The functions with hat are given by the previous expres- 
sions except for the propagator term which is removed by the amputation 
procedure. This is because the outgoing momentum must be on shelL 
The initial conditions are given by 

h^'^iPi) = eiiPi). A = ±1,0 
H{pi) = 1 

^ ~ I vx{-p^) if^,<0 

J/(„.) _ / ^A(Pi) iiEi>0 
^ ~ I M-Pi) if^i <o 

where the exphcit form of e^,u\,v\^u\^ v\ are given in the Ref. [20]. 

In order to actuahy solve these recursive equations it is convenient to 
use a binary representation of the momenta involved [10]. For a process 
with n external particles, to the momentum P^ defined as 

n 

P' = T.Pt (2-3) 

i=l 

a binary vector fh = (mi, . . . ,m„) can be assigned, where its components 
take the values or 1, in such a way that 

n 

pM = ^^.pM (2.4) 

i=l 

Moreover this binary vector can be uniquely represented by the integer 



m = 

i=l 



^2-im, (2.5) 
where 

1 < m < 2"-^ (2.6) 
Therefore all subamplitudes can be labeled accordingly, i. e. 

V'^(P) ^ V^(m), 



(2.7) 
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H{P) H{m). 

A very convenient ordering of integers in binary representation relies on the 
notion of level I, defined simply as 

n 

l = Y,mi. (2.8) 

i=l 

As it is easily seen all external momenta are of level 1, whereas the to- 
tal amplitude corresponds to the unique level n integer 2" — 1. This or- 
dering dictates the natural path of the computation; starting with level-1 
sub-amplitudes, we compute the level-2 ones using the Dyson-Schwinger 
equations and so on up to the level n which is the full amplitude. 

Contrary to original HELAC [20,21], the computational part consists of 
only one step, where couplings allowed by the Lagrangian defined by fusion 
rules are only explored. Subsequently, the helicity configurations are set up. 
There are two possibilities, either exact summation over all helicity config- 
urations is performed or Monte Carlo summation is applied. For example 
for a massive gauge boson the second option is achieved by introducing the 
polarization vector 

= e^'^4(p) + e-'^E't{p) + (2.9) 

where 4) G (0, 27r) is a random number. By integrating over cf) we can obtain 
the sum over helicities 

\=± 

The same idea can be applied to the helicity of (anti)fermions. 

Finally, the color factor is evaluated iteratively. Once again, we have two 
options. Either we proceed by computing all S""' x S"? color configurations, 
where the gluon is treated as a quark- antiquark pair and n^, Uq is the number 
of quarks and antiquarks respectively, or Monte Carlo summation is applied. 
Only a fraction of all possible S"' x S"' color configurations gives rise to 
a non zero amplitude. In the Monte Carlo approach for each event we 
randomly select a non vanishing color assignment for the external particles 
and evaluate the amplitTide. An overall multiplicative coefficient must be 
introduced to provide the correct normalization. The weight of the event 
is simply proportional to the |A^p multiplied by the number of non zero 
color configurations. Assuming that, on average, all color configurations 
contribute the same amount to the cross section this approach is numerically 
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more efficient than summing each event over aU colors, see Ref. [28] for 
further details. 

For the spinor wave functions as well as for the Dirac matrices, wc have 
chosen the 4-dimensional chiral representation which results in particularly 
simple expressions. All vertices in the unitary gauge have been included. 
Both the fixed width scheme (FWS) and the complex mass scheme (CMS) 
for unstable particles are implemented. 

The computational cost of HELAC grows like ~ S**, which essentially 
counts the steps used to solve the recursive equations. Obviously for large 
n there is a tremendous saving of computational time, compared to the n\ 
growth of the Feynman graph approach. 



3. Numerical Results 

As an example the algorithm has been used to compute total cross sec- 
tions for the production of weak vector bosons and the Higgs boson in 
association with jets. The following Standard Model input parameters have 
been used [29]: 

mw = 80.425 GcV, Tw = 2.124 GeV 

= 91.188 GcV, T^ = 2.495 GeV 

= 1.6637 X 10"^ GcV-5 

sin^ 9w = 'i- — rn^,/m%. 

The electromagnetic coupling is derived from the Fermi constant G^j^ ac- 
cording to 

V^G^mf^. sin^ Ow ,^ 

"cm = '^—^ • (3.1) 

vr 

All results are obtained with a fixed strong coupling constant as calculated 
at the mz scale 

as{ml) = 0.1187. (3.2) 



The mass of an intermediate and a heavy Higgs boson and associated 
Standard Model tree level widths are assumed to be: 

ruH = 130 GeV, Th = 0.005 GeV 

ruH = 200 GeV, Th = 1 GeV 

For the massive fermions the following masses have been applied: 
m„ = 4 MeV, = 8 MeV 
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Table 1. Results for the total cross section for the associated production of Higgs 
boson (130 GeV) with ati pair, CTbxact corresponds to summation over all possible 
color configurations, while CTmc corresponds to Monte Carlo summation. 



Process 


o-EXACT ± e (nb) 


(Tmc ± £ (nb) 


gg ^ ttH 


(0.2723 ± 0.0016) xlO~^ 


(0.2713 ± 0.0013) X 10--^ 


uu —>■ ttH 


(0.2758 ± 0.0017) X 10-* 


(0.2739 ± 0.0011) X 10--* 


dd ttH 


(0.1816 ± O.OOlljxlO-* 


(0.1811 ± 0.0007) X 10-^ 


cc ttH 


(0.8118 ± 0.0057) xlO"'' 


(0.8094 ± 0.0032) xlO"'' 


ss ttH 


(0.2203 ± 0.0014) xlO"'' 


(0.2191 ± 0.0008) X 10-^ 


bb ttH 


(0.2260 ± 0.0016) xlO"'^ 


(0.2262 ± 0.0009) X 10-'' 



Table 2. Results for the total cross section for the production of a heavy Higgs boson 

( 200 Ge V) via the vector boson fusion, (Texact corresponds to summation over all 
possible color configurations, while a^c corresponds to Monte Carlo summation. 



Process 


<7exact ± £ (nb) 


o-Mc ± £ (nb) 


uu uuH 


(0.140C ± 0.0029) xlO"-' 


(0.1361 ± 0.0020) X 10"^ 


uu —>■ ddH 


(0.6699 ± 0.0088) X 10-° 


(0.6596 ± 0.0081) X 10"^ 


ud — > udH 


(0.2280 ± 0.0043) xlO^'' 


(0.2222 ± 0.0021) X 10"* 


ud — > udH 


(0.1241 ± 0.0027) X 10"^ 


(0.1258 ± 0.0025) X 10"^ 


dd ddH 


(0.3404 ± 0.0046) xlO"^ 


(0.3477 ± 0.0032) xlO-'^ 


uu — uuH 


(0.5132 ± 0.0081) X 10"^ 


(0.5178 ± 0.0060) X 10-^ 



rris = 130 MeV, = 1.35 GeV 

mb = 4.4 GeV, 

mt = 174.3 GeV, Vt = 1.56 GeV 



The mixing of the quark generations is neglected. 

The CMS energy was chosen -y/s = 14 TeV. The following cuts were 
used to stay away from soft and collinear divergencies in the part of the 
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Fig. 3. Invariant mass distribution of the bb system in gg — > W+VI^ bbbb process. 



i 



Fig. 4. Transverse momentum distribution (left panel) and rapidity distribution 
(right panel) ofb-jet in gg W^W~bbbb process. 



phase space intergated over: 

PTi > 60 GeV, \yi\ < 2.5, Ai? > 1.0 (3.3) 

where 

Pr. = ^5CT^, ,. = lln(|±|l) (3.4) 

are the transverse momentum and rapidity of the i-jet respectively. Addi- 
tionaly Ai? is a radius of the cone of the jet defined as 

Ai?=^($,-$,.)2 + (y,-y,.)' 

where A$jj = $i — 

A#., = arccos { ^Mf.±ZA . 

V PTiPTj J 
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Fig. 5. Invariant mass distribution of the W'^W system in uu — > W^W dd 
process. 



Fig. 6. Transverse momentum distribution (left panel) and rapidity distribution 
(right panel) of d-jet in uu ~> W^W" dd process. 

There are several parameterizations for the parton structure functions, 
we used CTEC16 PDF's parametrization [30,31]. For the phase space genera- 
tion we used two different generators. The first one is PHEGAS [32], which 
automatically constructs mappings of all possible peaking structures of a 
given scattering process. Self-adaptive procedures like multi-channel opti- 
mization [33] is additionally applied exhibiting high efficiency. The Second 
one is a flat phase-space generator RAMBO [34]. 

In the Table ^ results for the total cross section for the associated 
production of the Higgs boson of mn = 130 GeV mass, with a tt pair 
are presented. The ttH production channel is one of the most promising 
reactions to study both the top quark and the Higgs boson at the LHC, in 
the second case especially for the bb decay channel of the Higgs boson [35,36]. 
As we can see from the Tabled the gg ttH process dominates due to the 
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Fig. 7. Transverse momentum distribution (left panel) and rapidity distribution 
(right panel) ofW in uu W^W~ dd process. 



enhanced gluon structure function. The final state of this channel consists 
of W bosons and four 6— jet, two from the decay of the top quarks and two 
from the decay of the Higgs boson. The main background process is gg — > 
W'^W~bbbb with the contributions from all intermediate states. In FigEl 
we present the invariant mass distribution for the bb system. Transverse 
momentum and rapidity distributions of 6-jet for the background process 
are shown in Fig. [l] 

A more powerful channel for higher Higgs boson masses is the vector 
boson fusion qq V*V* ^ qqH with H W^W decay. In the Table HI 
results for the total cross section for some production processes of a heavy 
Higgs boson of niH = 200 GeV mass via vector boson fusion are presented. 
The main background processes consist of qq — > W'^W^qq channels. As an 
example the distribution of the invariant mass of the W~^W~ system from 
the uu — > W~^W~dd process is presented in Fig|21 Transverse momentum 
and rapidity distributions of d-jet and W are also shown in Fig 1^1 and Fig|7| 



Summary 

An efficient tool for automatic computation of helicity amplitudes and 
cross sections for multi-particle final states in the Standard Model has been 
presented. Matrix elements and cross sections are calculated iteratively 
by Dyson- Schwinger equations. We are free from the task of computing all 
Feynman diagrams for a given process, which can become impossible for the 
large number of particles involved. The computationally expensive proce- 
dures of summing over color and helicity configurations have been replaced 
by Monte Carlo summation. At this stage, the code is able to compute scat- 
tering matrix elements and cross sections for hard processes. In the next 



12 



step we plan to make calculations of fully hadronic final states in pp and pp 
collisions possible. In particular, we wish to include the the emission of sec- 
ondary partons and translate the emerging partons into primordial hadrons 
e.g. by interfacing this package with codes like PYTHIA [37] or HERWIG [38]. 
This kind of multipurpose Monte Carlo generator will be of great interest 
in the study of the TeVatron, LHC or e+e^ Linear Collider data. 
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